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Abstract. We present a system of iV-coupled Lienard type nonlinear oscillators 
which is completely integrable and possesses explicit N time-independent and N time- 
dependent integrals. In a special case, it becomes maximally superintegrable and 
admits {2N — 1) time-independent integrals. The results are illustrated for the N = 2 
and arbitrary number cases. General explicit periodic (with frequency independent of 
amplitude) and quasiperiodic solutions as well as decaying type/frontlike solutions are 
presented, depending on the signs and magnitudes of the system parameters. Though 
the system is of a nonlinear damped type, our investigations show that it possesses a 
Hamiltonian structure and that under a contact transformation it is transformable to 
a system of uncoupled harmonic oscillators. 



PACS numbers: 02.30.Hq, 02.30.Ik, 05.45.-a 



Completely Integmhle N -Coupled Lienard Type Nonlinear Oscillator 2 
1. Introduction 

In a recent paper we have shown that the modified Emden type equation (MEE) with 
additional hnear forcing, 



X + ?,kxx + k X + \x = Q 



where over dot denotes diff^erentiation with respect to t and k and A are arbitrary 
parameters, exhibits certain unusual nonlinear dynamical properties [1]. Equation ([1]) is 
essentially of Lienard type. For a particular sign of the control parameter, namely A > 0, 
the frequency of oscillations of the nonlinear oscillator ([1]) is completely independent of 
the amplitude and remains the same as that of the linear harmonic oscillator, thereby 
showing that the amplitude dependence of frequency is not necessarily a fundamental 
property of nonlinear dynamical phenomena. In this case (A > 0) the system admits the 
explicit sinusoidal periodic solution 

Asm{ujt + 6) 
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1 - (^) AcosM + 5)' 
where A and 6 are arbitrary constants. In figure [Ih we depict the harmonic periodic 
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Figure 1. (Color online) Solution and phase space plots of equation llj for the case 
A < (a) Periodic oscillations (b) Phase space portrait 



oscillations of the MEE for three different initial conditions, showing the amplitude 
independence of the period or frequency. The phase space plot, figure [T)d, resembles that 
of the harmonic oscillator which again confirms that the system has a unique period of 
oscillations for A > 0. 

For A < 0, equation ([1]) admits the following form of solution [1], 




where Ji and I2 are constants. Depending on the initial condition, the solution is 
either of decaying type or of aperiodic frontlike type, see figure 2a. The time of decay 
or approach to asymptotic value is independent of the amplitude/initial value, which is 
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Figure 2. (Color online) (a) Decaying and frontlike solutions of ([T]) for the parametric 
choice A < 0, (b) Solution plot of H]) with A = 



once again an unusual feature for a nonlinear dynamical system [1]. Finally, for A = 
equation ([1]) is nothing but the MEE [3] which has the exact general solution [4] 

^^^^ " 2kt^ + hkt + h' 
where Ii and I2 are the two integrals of motion (figure 2b). 

A natural question which now arises is whether there exist higher dimensional 
coupled analogues of ([H) which are integrable and exhibit interesting oscillatory 
properties. In this paper we first report a system of two-coupled MEEs (with 
additional linear forcing) which is completely integrable, possesses two time-independent 
integrals and two time- dependent integrals and whose general solution can be obtained 
explicitly. Depending on the signs of the linear term, the system admits periodic 
(with amplitude independent frequency) or quasiperiodic solutions or bounded aperiodic 
solutions (decaying or frontlike). The results are then extended to A^-coupled MEEs and 
we prove the complete integrability of them also in the same way. From the nature of 
the explicit solutions we identify a suitable contact transformation which maps the 
coupled system onto a system of coupled canonical equations corresponding to a system 
of A^-uncoupled harmonic oscillators, thereby proving the Hamiltonian nature of the 
A^-coupled MEEs. We also prove that the system becomes a maximally superintegrable 
one [5] for any value of A^(> 1) when the coefficients of the linear force terms are equal. 

We organize our results as follows. In Section 2 we propose a two-coupled version 
of the MEE ([1]) and construct the integrals of motion through the recently proposed 
modified Prelle-Singer procedure (a brief description of this procedure is given in 
Appendix A). In Section 3, by appropriately choosing the coefficients of the linear forcing 
term, we construct different forms of general solutions including periodic, quasiperiodic 
and frontlike solutions and discuss the dynamics in each of the cases in some detail. In 
Section 4, we briefly analyse the symmetry and Painleve singularity structure properties 
of the proposed two-coupled version of the MEE. We also investigate the dynamics of 
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this equation under perturbation by numerical analysis. We consider a A^- coupled 
generalization of equation ([1]) and discuss the general solution/dynamics in Section 
5. In Section 6, we identify a Hamiltonian structure for the A^-coupled MEE by 
mapping it onto a system of iV-uncoupled harmonic oscillators through a contact type 
transformation which is obtained from the general solution of the coupled MEE. Finally 
in Section 7, we summarize our results. In Appendix A, we briefly discuss the generalized 
modified Prelle-Singer procedure which has been used to derive the results discussed in 
Section 2. 

2. Two-dimensional generalizations of the MEE 

It is of considerable interest to study the dynamics of higher dimensional versions of 
the MEE (II]). Recently Carinena and Ranada studied the uncoupled two-dimensional 
version of equation ([1]) [6], 

X + 'ikixx + k\x^ + Xix = 0, 

y + Shyy + kjy^ + \2y = 0, (5) 

where ki, k2, Xi and A2 are arbitrary parameters, and analyzed the geometrical 
properties and proved that the above system is superintegrable [6]. On the other hand, 
Ali et al. [7] have analysed a system of two-coupled differential equations, which is a 
complex version of ([1]) with x = y + iz and A = 0: 

y = -^yy - zz) - {y^ - 3yz'^), 

z = -3izy + yz)-i3yh-z'). (6) 

The above system of equations is shown to be linearizable by complex point 
transformation and from the solution of the linearized equation, a general solution for 
the equation ([6]) has been constructed [7]. 

In [2] we have pointed out that equation ([1]), and its generalization as well as 
the A^*'* order version can be transformed to linear differential equations through 
appropriate nonlocal transformations. In particular equation ([1]) under the nonlocal 
transformation U = x{t)e^o '=2;(-r)dr g^^g transformed to the linear harmonic oscillator 
equation U + XU = 0, where x and U are also related through the Riccati equation 
X = jjX — kx"^. Substituting the expressions for U and U and solving the resultant 
Riccati equation one can obtain the solution ([2]). Now searching for possible extensions 
to higher dimensions by considering a generalized nonlocal transformation of the form 
U = xe^o f(^i'^)M'r))dT ^ Y = yglo 9{x{r),y{T))dT ^ where U and V satisfy the uncoupled linear 
harmonic oscillator equations, U + XiU = and V + X2V = 0, we try to identify 
the forms f{x,y) and g{x,y) so that the transformation can be written as a system of 
parametrically driven Lotka-Volterra type equations or time-dependent coupled Riccati 
equations. In particular, with the choice f = g = kix + k2y, the transformation 
becomes a set of coupled time dependent Riccati equations, x = (^^x — kix"^ — k2xy^ , 

y = ( — kixy — k2y'^ ) . Consequently one obtains a system of two coupled MEEs 
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with additional linear forcing, 

X = —2{kix + k2y)x — {kix + k2y)x — {kix + k2y)'^x — \ix = 0i, 

y = -2{kix + k2y)y - {k^x + k2y)y - {k^x + k2yfy - Aay = (p2, (7) 

where fcj's and Aj's, i = 1, 2, are arbitrary parameters. When either one of the parameters 
ki or k2 is taken as zero, then one of the two equations in ([71) reduces to the MEE defined 
by ([1]) while the other reduces to a linear ordinary differential equation (ODE) in the 
other variable or vice versa. On the other hand when one of the variables {x or y) is 
zero, equation reduces to a MEE in the other variable. A characteristic feature of 
this form ([7]) is that it can be straightforwardly extended to higher dimensions as we 
see in the following sections besides admitting unusual nonlinear dynamical properties. 

To obtain the solutions of the above system of nonlinear ODEs one can solve the 
above coupled Riccati equations. However, to obtain the integrals of motion as well as 
the solutions we find it more convenient to solve ([7]) by the generalized modified Prelle- 
Singer (PS) procedure introduced recently [12]. We indicate this procedure applicable to 
(I7j) briefly in Appendix A. The resultant independent integrals of motion can be written 
as 

{x + {kix + k2y)xY + Aix^ 



h 



\{x + (fcix + k2y)x) + g(?/ + {kix + k2y)y) + 1 
{y + {kix + k2y)yf + Aa?/^ 



^^(i; + (fcix + k2y)x) + g(?/ + (fcix + k2y)y) + 1 

Ai > 
Ai <0, 



2 ' 



2 ' 



tan ^ 




- 


x+{k\x+k2y)x 




x+{k\x+k2y)x— 





x+{k-ix+k2y)x+^/\\-i\x 



tan ^ 


VMy 


- 


y+{k\x+k2y)y 




[y+{kix+k2y)x- 





A2 > 
A2 < 0. 



(8) 
(9) 

(10) 
(11) 



y+(k-^x+k2y)y+^J\\2\y 

We note here that the forms of the time-dependent integrals of motion depend 
upon the signs of the parameters Aj, i = 1,2. Using the above integrals one can obtain 
periodic and aperiodic but bounded solutions as per (i) Ai, A2 > (ii) Ai,A2 < and 
(iii) Ai < 0, A2 > (or vice versa). The case Ai = A2 = is dealt with separately below 
in Section 3.3. In the following we discuss the nature of the solutions. 



3. The Dynamics 

3.1. Periodic and quasi periodic oscillations (Ai, A2 > 0) 

By restricting Ai, A2 > in the integrals (l8|)- (fTT|) and solving them algebraically we 
obtain the following general solution, 

A sin(co'it + 61) 



x{t) 



1 _ ^ cOs(^it + COs{u2t + 62) 
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y{t) 



B sm{ijj2t + ^2) 



^ COs(cU2t + 62 



Aki Bk2 



L02 



< 1, 



(12) 



Aj, j = 1,2, A = vA/t^i, B = V^/u;2, 5i = I3, 62 = h- Two types 



where Uj — y 

of oscillatory motion can arise depending on whether the ratio uji/uj2 is rational or 
irrational leading to m : n periodic or quasiperiodic motion, respectively. One may 
note that the frequency of oscillations is again independent of the amplitude in the 
present two-coupled generalization also. The conservative nature of the above oscillatory 
solution can be seen from the phase space plot. 
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Figure 3. (a) Quasi-periodic oscillations with uji — 1 and lj2 = V2 (i) phase space 
plot (ii) Configuration space plot and (iii) Poincare SOS (b) periodic oscillations with 
ciJi = 1 and UJ2 — 5 (i) phase space plot (ii) Configuration space plot and (iii) Poincare 
SOS, see equation ([T! 



To visualize the dynamics we plot the solutions x(t) and y{t) both in phase space 
and in configuration space for two different sets of frequencies. First we consider the 
case in which the ratio of frequencies is an irrational number. For illustration, we 
take uji = uj2/\/2 = 1. The {x,x) phase space plot is depicted in figure 3a(i). The 
configuration space (x, y) plot is given in figure 3a(ii). We also confirm the quasi-periodic 
nature of the oscillations by plotting the Poincare surface of section (SOS) [13] in figure 
3a(iii). The solution is clearly not a closed orbit and is quasiperiodic or almost periodic 
in the sense that trajectory returns arbitrarily close to its starting point infinitely often, 
which is confirmed by a closed curve in the Poincare SOS. The system parameters are 
chosen as fci = — 1, /c2 = 1 and the values of the arbitrary constants are taken as A = — 1 
and B = -0.4. 

Next, when the ratio of frequencies is a rational number one has periodic solutions. 
In figure 3b we show a 1:5 periodic solution by choosing ui = 1 and 002 = 5. 
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3.2. Aperiodic solutions (Ai,A2 < 0) 

We next consider the case Ai < 0, A2 < 0. In this case one obtains aperiodic but bounded 
fronthke or decaying type solutions. To see this, we use the two time-independent 
integrals (hjh) and time-dependent integrals {I3, h) vide equations (l8l).(!9!). (fT0|) and 
( ITTl) . and obtain the following general solution for equation ([7]) when Ai, A2 < 0, 

x{t) 



lAilt^, 



|Ai|t 



+ 336- 



|Ai|t 



) + fc2Aip2eVl-^2|t + g^e 



\^2\t 



AivWpae 



/ciA2pieVl^il* + 'j3e 
where 3i = \/h/ {\/hXiX2), ^2 



+ fc2Ai(D^2e 



\X2\t 



(13) 



y72/(y74AiA2), 33 = v/7^/(AiA2), 3^ = 
y/hh/ (A1A2). Here also we observe that the solution is of the same form as in the one 
dimensional case. The decaying type solution is plotted in figure 4a and the frontlike 
solution is shown in figure 4b. 




t 6 
(a) 



{!>) 




Figure 4. (a) Decaying type solution of ^ for Ai, A2 < (b) Frontlike solution of 
([7]) for Ai, A2 < (c) Decaying type solution of ^ for Ai = A2 = 0. 



3.3. Decaying type solution for Ai = A2 = 

Restricting the values of the parameters to Ai = A2 = in ([7]) one can get the two- 
dimensional generalization of the MEE ([1]) . In this case we find that the system ([7]) 
admits the following integrals of motion, namely, 

(x + {kix + k2y)x) r _ 4 , ^ 
-'1 — . , /, ~~; — ^ — ) I2 — —t + 



y + {kix + k2y)y ' x + {kix + k2y)x^ 

/3 = _t + y h = -+ 1 - tjk^x + k2y) ^^^^ 

y + {kix + k2y)y' 2 {kix + k2y + {kix + k2yYy 
Then the general solution for equation ([7]) (with Ai = A2 = 0) can be written in 
the form 

2/l(/2 + t) 



x{t) 



kJi{2h + (2/2 + t)t) + A;2(2/4 + (2/3 + t)t) ' 



^ 2(/3+t) 

kih{2h + {2l2 + t)t) + k2{2h + {2h + t)t)' ^ ^ 
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Again the solution turns out to be a rational function in t in which the denominator 
is a quadratic function of t as in the one dimensional case. Choosing Ji, /2, 1^ and 1^ 
suitably, one can obtain decaying type solutions as shown in figure 4c. 



3.4- Cases with mixed signs of parameters Ai and A2 

In the mixed case, for example with Ai > 0, A2 < 0, the solution can be derived from 
(USD as 

—uiA sin(ci;it + 5i) 



x{t) 



yit) 



(j^eV l^^l* + J4e-V l^^l*) + Aki cos(cJit + 5i) - 1 



Ai, 



(16) 



2[^(J2eV l^^l* + J4e-v l^^l*) + Aki cos{cu^t + 5i) - 1] 
where A = \fTxlio\. 

The motion now turns out to be a mixed oscillatory-bounded frontlike one. We 
depict the solution (fT6l) in figure 5 in which we have fixed /ci = /c2 = 1, = 2, and 
A2 = -1. 





Figure 5. Solution plots of equation ([7|) for the mixed case Ai > 0, A2 < (a) a;(i), 
(b) 2/(0 



3.5. A superintegrable case 

Finally we investigate the dynamics in the limit Ai = A2 = A 7^ 0. In this case one is 
able to find an additional time-independent integral of motion given by 

{kix + k2y + {kix + k2yY + A)^ 
As a result one has three time-independent integrals, /i, /2, and /s, besides one time- 
dependent integral (any one of the two time-dependent integrals J3 or J4), in a two 
degrees of freedom system which in turn confirms that the system under consideration 
is a superintegrable one. Obviously it is also maximally superintegrable [5]. 
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4. Symmetry and Singularity structure analysis 

In order to understand why the system ([71) is integrable and whether perturbations 
of them lead to chaos, we analyse the symmetry properties, in particular the point 
symmetries and Painleve singularity structure associated with equation ^ and 
investigate numerically the perturbed equation. 



Lie point symmetry analysis 

Considering the invariance of equation (JTj) under a one parameter continuous Lie point 
symmetry group [8], a; ^ X = x + erjiit, x, y) + O(e^), y ^ Y = y + eri2{t, x, y) + O(e^), 
t ^ T = x + e^{t,x,y) +0(e^), e << 1, and performing the Lie point symmetry analysis 
on the equation ((71) using the package MULIE [9], we find that the system admits only 
the following single Lie point symmetry vector, 

corresponding to the invariance of ([7]) under time translation for Ai 7^ A2. For the 
specific choice Ai = A2, which we have earlier proved to be superintegrable, we obtain 
the following two additional point symmetries, 

ki d d 

k2 oy ox 

^ ko d d , ^ 

ki ox oy 

So we conclude at this stage that since the system ([7]) is completely integrable and 
does not possess enough number of Lie point symmetries, it has to admit more general 
symmetries, namely nonlocal and contact symmetries. The study requires separate 
analysis and we do not pursue them here. Note that on the other hand, the one 
dimensional MEE ([1]) admits eight Lie point symmetries and is linearizable through 
point transformation which is not the case for the coupled MEE ([7]). 



4-2. Painleve singularity structure analysis 

We now perform the standard Painleve singularity structure analysis [10, 11] to the 
equation ([7j). Looking for the leading order behaviour of the Laurent series solution in 
the neighbourhood of a movable singular point to, we substitute x = qot^ and y = boT'^, 
T = (t — to) ^ 0, and obtain 

aopip - l)rP^2 ^ ^3^2^3p ^ SalkipT^P-^ + 2aobok2pT'+''-^ + aobok^qr^^''-^ 

+ 2albokik2T^P^'' + aoblklr^+^^ = 0, (20) 

boqiq - l)r^-' + fej^fc^r^^ + aobok,prP+'^"' + 2aobok^qTP+'^-^ + alboklr^^+'i 

+ Sfo^fesgr^-?"! + 2ao&o^iA;2rP+'' = 0. (21) 
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Comparing the exponents of r, we find p = —1 and g = — 1. Substituting tliis and 
simplifying we get 

(2ao — SoqA;! + a^k\ — 80060^2 + 2ao6o^i^2 + 0.0^0^2)''""^ = O5 
(26o - 3ao&oA;i + alb^kl - 3blk2 + 2aoblkik2 + b^k^T^^ = 0. 

Solving the above system of equation we find two possibilities for the leading order 
coefficients as = (1 — &o^2)/^i and ag = (2 — &o^2)/^i5 while 60 is arbitrary. For 
the choice ao = (1 — ^0^2)/^! we identify that the resonances (that is powers at which 
arbitrary constants can enter) occur atr = — l,r = 0,r = l and r = 1. By proceeding 
with the full Laurent series one can show that in addition to to and 60 being arbitrary 
(corresponding to r = —1 and r = 0), Oi and bi are also arbitrary corresponding to 
r = 1, 1, while all higher order coefficients in the Laurent series can be determined in 
terms of the earlier ones. Similarly we find that the Laurent series corresponding to the 
second value of Oq = (2 — 60^2)/^! also does not admit any movable critical singular 
point. We find that the equation ([71) passes the Painleve test as expected. 

4-3. Perturbed system : Numerical analysis 

In order to understand the dynamics of the system in the neighbourhood of the integrable 
parametric regime, we perturb the system leading to the form 

X + 2{kix + k2y)x + {kix + k2ii)x + {kix + k2yYx + Aix + pixy = 0, 
y + 2{kix + k2y)y + {kix + k2y)y + {hx + k2yfy + A2I/ + P2xy = 0, (22) 

where pi, p2 are the strength of the perturbation. We now numerically solve equation 
( I22I) for three different values of pi and p2 using the fourth order Runge-Kutta method. 
Of these three parametric choices, the first one (pi = P2 = 0) corresponds to the 
completely integrable equation (I7j) and the other two are the perturbed form of ([7]), 
namely, pi = p2 = 0.5 and pi = P2 = 0.7. In figures 6a we show the phase space plots of 
equation fl22|) for these three different parametric choices. The corresponding Poincare 
surface of sections are plotted in figures 6b by identifying the peaks of y{t) from the 
time series data. We find from the Poincare surface of section of the torus that for 
the integrable parametric choice a closed curve is obtained while as the strength of the 
perturbation is increased the curve starts to break up and ends up into scattered points, 
confirming the onset of chaos. 

5. A^-coupled MEE 

Proceeding further, we find that equation ([T]) can also be generalized to arbitrary 
number, A^, of coupled oscillators. In order to generalize the results of two-coupled 
MEE ([7]) to A^-coupled MEE, we have first investigated the three-coupled MEEs. Then 
from results of two and three-coupled MEE, we have generalized the results to A^-coupled 
MEE. For brevity, we are not presenting here the results of the = 3 case. 
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Figure 6. Figures (a) (i), (ii) and (iii) describe phase space trajectories of equation 
with Ai = 1 and A2 = 2 for three values of pi, pi corresponding to integrable 
case (i) p\ = P2 = ^ and perturbed cases : (ii) p\ — pn — 0.5 and (iii) p\ — p2 = 0.7. 
Figures (b) (i) , (ii) and (iii) describe the Poincare surface of sections of the above phase 
space plots. 



The A^-coupled modified Emden type equation is given as 

N N N 

Xi + 2(^^ kjXj)xi + (^^ kjXj)xi + kjXjYxi + AjXj = 0, i = 1, 2, A^. (23) 
i=i i=i j=i 

By generalizing the results of two and three-coupled cases one can now inductively 
construct the integrals of motion for the equation (l23|l which turn out to be 

(xj + '^j—i{kjXj)xiY' + 



'2i 



+ 1 



2 ' 



tan 



\jXi 



Xi -\- ^y^j j—] i,kjX j)Xi 

In other words one has N time-independent integrals and time-dependent integrals. 
From these integrals one can derive the general solution of equation (123|1 for the 
parametric choice Ai, A2, Atv > as 

N 



Xit, Xi> 0, i 



l,2,...,iV. 



(24) 



(25) 



1 - Ef=i ^ cos(^j-t + 5j) ' 



UJj 



< 1, 



(26) 



with uji = \/Xi. 

In the second case, Aj < 0, i = 1,2, ...,N,to derive the general solution we construct 
the time-dependent integrals in the form 



Xi -\- ^^j—ii^kjXj^Xi \Xi\xi 



l3i 



Xi ~\~ ^ -^ j'=1 i^i^ i^^i "v I 



1,2, ...,iV. 



(27) 
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Using the above integrals one can construct the general solution of fl23l) in the form 



Xi{t) = iis=i ^Vi ^\K 3^ J ^ z = 1, 2, N, {21 



where j and s 7^ Jii = ^/hil {Vhi]\n=i ^n), ^si = ^hihil X\n=\ ^n- 

In the third case we consider the mixed sign case of the parameters Aj, i = 1, 2, A^. 
Let us consider that A^, r = 1,2, ...,l, are positive constants and Xj, j = N — I, N, 
are negative constants. The solution for this case can be written as 

Xr = —UrAr sin(a;r^ + Sr)h~'^, 
I N 

^^• = n^' n A„y^pi,eV^* + J3,eV^*)(2/i)-\ (29) 

i=l m=N-l 

where 

/ N I N 

h = J2krArCOs{uJrt + 6r)+ ^™ H "^-^ 11 Agpi^cV^* + Jg^e-V^*) - 1, 

r=l m=N-l r=l q=N-l 

m ^ q and m ^ j. 

In the fourth case we have the A^-dimensional generalization of "classical MEE". 
It is evident from the lower dimensional cases that the integrals of motion includes (i) 
(A^ — 1) time-independent ones (ii) time-dependent integrals which are linear in 't' 
and (iii) one time- dependent integral which is quadratic in 't', that is 

/n = ,^^^-^%^, . = 1,2,...,A^-1, (30) 

I,^ = -t + - t = l,2,...,N, (31) 

2 (E,=i kjXj + {kjXjy) 



From these integrals, again one can deduce the general solution of (1231) for Aj = in the 
form 

2Iu{l2i + t) 

Y.''k,h,{2h + {2h,+t)t) 



x^{t) = _^ , ^ 1 z = l,2,...,Ar-l, (33) 



Mt) = -1, ^Shulll , (34) 

E7=iVy(2/3 + (2/2,+t)t) 

where Ii^ = 1. 

In the fifth case, one can consider the parameters Aj's are all equal but nonzero, 
that is Ai = A2 = ... = Atv = A 7^ 0. In this case, one can construct (A^ — 1) additional 
time-independent integrals by eliminating the variable 't' in equation (1251) . 



'3. = n'":-"] ^ = 1, 2, ...A^ - 1. (35) 



TijikjXj + {kjXjY + A) 
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Again the existence of (2A^ — 1) time-independent integrals of motion (vide equations 
flM|) and fl35|) ) confirms that the system under consideration, namely fl23|) with Aj = 0, 
i = 1,2, N, is a maximally superintegrable one. 



6. Connection to uncoupled harmonic oscillators 

The solution of the coupled MEE equation (171) given by equation (JT2l) and (fT3i) can be 
rewritten as 

, = Z ^ (36) 

Ujf U)f 

where U = A sin(co'it+5i) and V = B siia{uj2t+62). Here U and V can also be interpreted 
as the solutions of the following uncoupled harmonic oscillator equations, 

t) + ujIU = 0, V + ulV = 0. (37) 

Equation fl36l) gives a transformation connecting equation ([7]) and the harmonic 
oscillator equations ( !37l) . In order to invert this transformation, we need U and V for 
which we differentiate equation fl5Bl) once with respect to time and replace U, V with 
—ufU and — 1^2^ respectively. Thus we get the following equations 

tin - ^ - ^) - uikiU + k2V) vii-^-^)- vikiU + k2V) 

Solving the relations (1361) and (I38p one obtains the inverse transformation of (I36p and 
(|38|) as 



U = ..... , , ..... ,1... . , . .. . (39) 



^(i + {kix + k2y)x) + %y + [kix + k2y)y) + 1) 



"2 



y = 7177-^^ . , . . .1.. . . , . . . (40) 



■^(x + (fcix + fc2y)a;) + z^{y + {kix + k2y)y) + 1) ^ 



(i + kix"^ + /c2a;?/) 
r(i; + {kix + fc2Z/)a;) + ;^(2/ + (^la; + k2y)y) + 1) 



rj ^ -r -r 1^2-^y) 



y ^ jy + hxy + k2y^) ,^2) 

(^(i + (kix + k2y)x) + ^(2/ + (fcix + k2y)y) + 1) ' 

This is indeed a contact transformation between the old and new variables, which 
can also be interpreted as a linearizing transformation to equation ([7]). The form of the 
Hamiltonian for the system of two uncoupled harmonic oscillators ( 137]) obviously is 

H = ]^ [Pi + P| + Aif/' + \2V^] , (43) 

where Pi = U, P2 = V, Ai = and A2 = Substituting for Pi, P2, U, and ^ from 
([39])- (02]) we get, 

^ ^ (x + fcix^ + ^2X1/)^ + {y + kixy + k2y^)^ + Aix^ + Asy^ 

2[f^(x + (fcix + A;22/)x) + 1(1/ + (A;ix + My) + 1]' 
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Here we note that H = + 12)-, where Ji and I2 are the time independent integrals 
of motion of equation ([7]) (vide ([8]) and ([9])). 

Further, U and V satisfy the canonical equations 

SH ^ . 3H ... ^^gj 



U 

V 



dPi 
dH 



Pi, Pi 



P., P, 



dU 
dH 



-Aif/, 
-X2V. 



(46) 



dP2 dV 
The results obviously confirm the existence of a conservative Hamiltonian for equation 

The above results can also be extended to dimensions by using the contact 
transformation derivable from the solution ( l26i) . 



Xi -\- (^ ^j— 1 



Pi = 



+ 1 



N 



+ 1 



(47) 



z = 1,...,A^. Consequently the Hamiltonian for equation fl23|) can be rewritten as a 
system of uncoupled harmonic oscillators specified by 

N 

(48) 



1 ^ 



i=l 



In terms of the original coordinates this becomes 

{xi + i^^j—i kjXj)xi)'^ + XiX^ 



H 



2^j=i 



kj / ■ 



1 2 ' 



+ 1 



(49) 



which can be associated with the integrals of motion flM|) as if = | X]j=i -^li- The 
canonical equation of motion are given as 



Ui 



dH 



P 



dH 



1,2,..., AT, 



(50) 



dP,' ^' dU 

thus confirming the Hamiltonian nature of the system ( l23l) . However, we have not yet 
succeeded to obtain an explicit Lagrangian form, and so canonically conjugate momenta, 
to re-express ( 144|) or ( l49l) in terms of canonical coordinates. This is being pursued at 
present. 



7. Conclusion 



In this paper we have presented a system of completely integrable A^-coupled Lienard 
type (modified Emden type) nonlinear oscillators. The system admits in general A^ time- 
independent and A^ time-dependent integrals whose explicit forms can also be found. For 
special parametric choices, the system also becomes maximally superintegrable. Using 
these integrals general solution of periodic, quasiperiodic, frontlike and decaying type or 
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oscillatory type are obtained depending on the signs and magnitudes of the linear forcing 
terms. We have also pointed out that the system possesses a nonstandard Hamiltonian 
structure and is transformable to a system of uncoupled harmonic oscillators. Further 
analysis of the Hamiltonian structure can be expected to yield interesting information 
on the nonstandard Hamiltonian structure of coupled nonlinear oscillators of dissipative 
type. It is also of interest to investigate whether there exist other couplings of MEEs and 
its generalizations which are also integrable : For example, one can show [2] that under 
the general transformation U{t,x) = x"e'^o ■/'(^(*'))'^*', f{x{t)) = (3{t)x"^ + 7(t), the one 
dimensional generalized MEE, x + {n — l)^ + ^x'^'^^^ + bi{t, x)x + b2{t)x"^^^ +b3{t)x = 0, 
where bi{t,x) = ^ (2^7 + nA + (m + 2n)/3x"), = + 27/3 + A/?) , = 

^(7 + 7^ + A7), can be reduced to the linear harmonic oscillator equation. One can 
expect higher dimensional generalization of such a transformation can give rise to more 
general higher dimensional integrable equations. These questions are being pursued 
currently. 
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Appendix A. Generalized modified Prelle-Singer procedure 

To solve the system of two-coupled second order nonlinear ODEs we apply the 
generalized modified Prelle-Singer (PS) approach introduced recently [12]. Let the 
system 

X = —2{kix + k2y)x — {kix + k2y)x — {kix + k2yYx — Aio; = 0i, 
y = -2{kix + k2y)y - {kix + k2y)y - {kix + k2yfy - Aay = 02, 

admits a first integral of the form /(t, x, i, y) = C with C constant on the solutions 
so that the total differential gives 

dl = Itdt + I^dx + lydy + Ij.dx + lydy = 0. (A.l) 

Equation ([7]) can be rewritten as the equivalent 1-forms 

4)idt - dx = 0, (f)2dt - dy = 0. (A. 2) 

Adding null terms Si{t,x,y,x,y)xdt — si(t,x,y,x,y)dx and S2(t,x,y,x,y)ydt — 
S2(t,x,y,x,y)dy with the first equation in ( 1A.2I) . and ui(t,x,y,x,y)xdt — ui(t,x,y,x, 
y)dx and U2it,x,y,x,y)ydt — U2it,x,y,x,y)dy with the second equation in (1A.2I) . 
respectively, we obtain that, on the solutions, the 1-forms 

(01 + six + S2y)dt — sidx — S2dy — dx = 0, (A. 3) 

(02 + uix + U2y)dt — uidx — U2dy — dy = 0. (A. 4) 
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Hence, on the solutions, the 1-forms (lA.ip and (lA.3p - (lA.4p must be proportional. 
Multiplying flA.3|l by the function R{t, x, y, x, y) and flA.4|) by the function 
K{t,x,y,x,y), which act as the integrating factors for ( ]A.3I) and ( ]A.4p . respectively, 
we have on the solutions that 

dl = i?(0i + Sx)dt + if (02 + Uy)dt - RSdx - KUdy - Rdx - Kdy = 0, (A.5) 

where S = {Rsi + Kui)/R and U = {Rs2 + KU2)/ K. Comparing equations (1A.5I) and 
flA.l|) we have, on the solutions, the relations 

It = Rih + Sx) + if (02 + Uy), I, = ~RS, ly = -KU, h = -R, k = -K. (A.6) 

The compatibility conditions between the different equations in ( 1A.6P provide us the ten 
relations 



D[S] 

D[U] 
D[R] 
D[K] 

SRy 
Ry 



>lx 



')2y 



K 

r' 

R 



K 

R' 
R 



>2x + sS<p2x + 5*01^ + S 



..^ly + —U(f)ly + U(f)2y + U' 



-{R(Pi, + K^2± + RS), 

-{K^2y + R^ly + KU), 

-RSy + UK, + KU,, R, 
UK., + KUi, 



UKy + KUy, 



Ry 



- SR± + RSx, 
a 



(A.7) 

(A.8) 
(A.9) 
(A.IO) 
(A.ll) 
(A.12) 
(A.13) 



Here the total differential operator, D, is defined by i^ = ^ + x-^ + y-^ + '/'i^ + 
Integrating equations flA.Op . we obtain the integral of motion. 



/ = n + r2 + rs + r4 - 



K + 4^(ri + r2 + r3 + r4^ 
dy 



dy, 



(A. 14) 



where 
ri = 



J (i?(0i + Sx) + if(02 + t/y)) dt, r2 = - J (i?5 



d 

dx 



(ri) jdx 



r3 = - 



KU + ^(ri + r2] 
dy 



dy, 



r4 = - 



d ( 

it + — (^ri + r2 + ra 



dx. 



Solving equations flA.7p - flA.13p one can obtain S,U,R and K. Substituting these 
forms into flA.14p and evaluating the resulting integrals one can get the associated 
integrals of motion. Once sufficient number of integrals of motion are found (four in the 
present problem) then the general solution can be derived from these integrals by just 
algebraic manipulations. One can refer to [12] for details of the method of solving the 
determining equations flA.7P -f lArT3l) . 
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